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ABSTRACT 

The inverse Compton (IC) interpretation of the excess EUV emission, 
that was recently reported from several clusters of galaxies, suggests that the 
amount of relativistic electrons in the intracluster medium is highly significant, 
W e > 10 61 erg. Considering Coma as the prototype galaxy cluster of nonthermal 
radiation, with synchrotron and IC fluxes measured in the radio and EUV 
regions, and possibly also in the hard X-ray region, we discuss implications of 
the inverse Compton origin of the EUV fluxes in the case of low intracluster 
magnetic fields of order 0.1 /xG as required for the IC interpretation of the 
observed excess hard X-ray flux, and in the case of high fields of order 1 /iG 
as suggested by Faraday rotation measurements. Although for such high 
intracluster fields the excess hard X-ray fluxes will require an explanation other 
than by the IC effect, we show that the excess EUV flux can be explained by 
the IC emission of a 'relic' population of electrons driven into the incipient 
intracluster medium at the epoch of starburst activity by galactic winds, 
and later on reenergized by adiabatic compression and/or large-scale shocks 
transmitted through the cluster as the consequence of more recent merger 
events. Radiative cooling will naturally produce a sharp cutoff in the spectrum 
of these relic electron population, which is required, in the case of /i-Gauss 
fields, in order to avoid a contradiction with the observed radio fluxes. For high 
magnetic fields B > 1 /iG the interpretation of the radio fluxes of Coma requires 
a second population of electrons injected recently. They can be explained as 
secondaries produced by a population of relativistic protons. We calculate the 
fluxes of 7-rays to be expected in both the low and high magnetic field scenarios, 
and discuss possibilities to distinguish between these two principal options by 
future 7-ray observations. 
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1. Introduction 

One of the significant characteristics of galaxy clusters is their nonthermal particle 
content. A clear signature of the presence of relativistic electrons, and by implication also 
of relativistic protons, in the intracluster medium (hereafter ICM) is the detection of a 
diffuse cluster-scale radio emission from a number of these objects ( fcee e.g. Giovannini et al 



1999| , and references therein). The radio fluxes give an important information on the energy 
spectra of the electrons in the ICM, but being dependent on the ICM magnetic fields, 
they do not yet allow unambiguous estimates of the overall energy content of relativistic 
electrons. More accurate estimates can be derived if one interprets the 'excess' diffuse 
EUV emission reported from some clusters ( [Lieu et al. 1996| , Mittaz et al. 1997| , |Bowyei| 



|fe Berghofer 1998] ; but see also [Arabadjis fc Bregman 1999|) in terms of inverse Compton 



(IC) upscattering of 2.7 K microwave background radiation by electrons with energies of 
a few 100 MeV ( |e.g. Hwang"T997i |Sarazin fc Lieu 1998| , |Enfflin fc Biermann 1998j) . This 



typically requires an overall energy in these electrons of about 5 x 10 61 erg. Depending on 
the spectra of relativistic electrons at lower energies, and on the ratio of the relativistic 
protons to electrons, the total energy in cosmic rays (CRs) might be then as high as 
W C r ~ (10 62 - 10 63 )erg. 

Thus, CRs may contain a significant fraction of the overall kinetic energy in the 
ICM and could be dynamically important. Because the energy loss time of the dominant 
hadronic CR particles in the ICM typically exceeds the Hubble time, and since their 
confinement time generally exceeds the cluster lifetime ( |Volk et al. 1996| , [Berezinsky et al. 



1997| , |Colafrancesco fc Blasi 1998|) , the nonthermal component of galaxy clusters should 



contain valuable information about the evolution of these largest gravitationally bound 
systems of the Universe. 

In this paper we study in detail the broad band nonthermal emission of the Coma 
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cluster, considered as a prototype, with nonthermal radiation fluxes reported both in the 
radio and the EUV regions. In the Coma cluster the 'excess' (over a single-temperature 
plasma emission) flux is firmly established also in the region of hard X-rays by |Fusco- 
[Femiano et al.| (1999) and [Rephaeli et al.| (1999). The interpretation by these authors of 
this hard X-ray emission in terms of IC radiation of electrons which are also responsible for 
the radio fluxes at low frequencies, suggests that the ICM magnetic field in Coma is about 
0.1 — 0.2 /zG. In section 2 below we summarize the results following from the hypothesis 
that, along with the excess EUV flux, the excess hard X-ray flux also has an IC origin, and 
calculate the fluxes of IC radiation to be expected in the 7-ray region. 

However, it is quite possible that the ICM magnetic field in Coma is at a level of 
> 1 /uG, as deduced from Faraday rotation and depolarization measurements on background 
radio sources ( [Kim et al. 199ti| , [Feretti et al. 1995| ). Although an IC origin of the excess 
hard X-ray flux in that case is excluded, we show in section 3 that an IC origin of the excess 
EUV flux is still possible. We demonstrate that a rather unusual energy distribution of the 
relativistic electrons, with a sharp cutoff above a few 100 MeV (that is needed to comply 
with the data both in EUV and the radio region), can be formed during the cosmological 
evolution of the Coma cluster. In particular, these electrons could be the relics of the 
early starburst activity in the Coma galaxies which fed the incipient ICM with relativistic 
particles and magnetic fields via their galactic winds ( |Volk fc Atoyan 1999| ). In section 4 
we calculate the 7-ray fluxes to be expected in the high magnetic field case and discuss 
whether it would be possible to distinguish between high and low magnetic field scenarios 
by future 7-ray observations. 
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2. IC interpretation of the excess hard X-rays 
Radio fluxes from the Coma cluster are detected from 10.3 MHz to 2.7 GHz ( [Bridle 



Purton 196S| , Kenning 1989| , |Kim et al. 1990| , |Giovanuini et al. 1993|) , and at 4.85 GHz an 



upper flux limit was reported by Kim et al. (1990). The differential energy fluxes J u oc v~ ar 
between 30.9 MHz and 1.4 GHz are well fitted with a power-law index a r = 1.16. The data 
at 2.7 and 4.85 GHz fall below this extrapolation, so that the 'best fit' single power-law 
index would then correspond to a r = 1.36 flSchlickeiser et al. 1987| ). Note however that 
the reported sharp decline of the radio fluxes above 1.4 GHz can also be explained as an 
instrumental effect (Deiss 1997). 



The diffuse excess (over the thermal) EUV radiation between 70 eV and ~ 200 eV has 
a mean power- law index a euv ~ 0.75 (Lieu et al. 1999). The power- law index of the diffuse 
excess X-rays is rather uncertain and lies most probably within the range of a x ~ (1 — 2) (see 
Fusco-Femiano et al. 1999; Rephaeli et al. 1999), which includes the values of the spectral 
indices a T of the measured radio fluxes. The latter circumstance becames important for an 
IC interpretation of the excess X-ray flux of Coma because this radiation should be then 
produced by the same population of electrons which are responsible for the synchrotron 
emission at radio frequencies. 

The characteristic frequency of synchrotron radiation of electrons with Lorentz factor 



7 is v ~ _B At G7 2 Hz where B^q = B/ljiG Qe.g. Ginzburg 1979|) . Thus the production of 



synchrotron radiation in the region from 10 MHz to 1.4 GHz requires that the energy 
distribution of the electrons N(j) cx r y~ a ° extends with an index a e = 1 + 2a T rj (3.2 — 3.6) 
from 7 r ,i ~ 3 x 10 3 B~g 5 to 7 r , 2 ~ 4 x 10 4 B~q 5 . 

For the IC process the mean energy of photons produced in the 2.7 K background 
radiation field with the mean photon energy eo — 3 kT is e ~ (4/3)eo7 2 — 8.7 x 10 _4 7 2 eV. 
For the production of the X-ray photons from 20 keV to 80 keV one therefore needs electrons 
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with energies between 7 Xjl ~ 



5 x 10 3 and 7 Xi2 ^ 10 4 . Because the latter actually corresponds 



to the lower edge j Tj i of the electrons responsible for the observed radio emission (in 
magnetic fields of order 0.1 /xG), the spectral index of the X-ray emitting electrons should 
be close to a c as derived from radio data. 

The IC origin of Coma's EUV spectrum from 70 eV to 200 eV implies a spectral index 
a e ~ 2.5 for the electrons in a rather narrow energy region around 7 CUV ~ 300 (Sarazin & 
Lieu 1998). A gap by a factor of about 10 between 7 euv and 7 rjl allows then a scenario 
with a single population of power-law electrons where radiative losses induce a break in a c 
somewhere in that gap. Thus assuming a source function of electrons Q(j) oc j~ ain i in the 
ICM with a in j = 2.3, the IC emission of those electrons at low energies will give a euv = 0.65, 
which is acceptably close to the value 0.75. After the break, the electron distribution 
steepens to a e = 1 + a^j = 3.3 required for interpretation of the radio fluxes. 

Such a break in the energy spectrum can be produced by radiative losses of electrons if 
they are injected into the ICM on time scales At in j > 10 9 yr. For the radiative energy loss 
rate due to IC (in the Thompson limit) and synchrotron processes (e.g. Ginzburg 1979) 



where we take into account the cosmological evolution of the energy density of the cosmic 
microwave background radiation (MBR), u> m b r oc (1 + z) 4 at the epoch z. For magnetic 
fields B < 3/iG the radiative losses are dominated by the IC process. Using Eq. (2), from 
the condition t ra d(7) — At in j we find that the radiative break would appear in the 'proper' 
place (7 br ~ 10 3 ) if relativistic electrons were injected into the ICM recently, during times 
£inj ~ (1 — 3) x 10 9 yr (this corresponds to z- m j < 0.2). Note that the relativistic electrons 




(1) 



the characteristic energy loss time is 



t rad » 2.4 x 10 12 [(1 + zf + O.ISJg]" 1 7" 1 yr , 



(2) 
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need to be produced/accelerated in the ICM continuously during that time period, because 
in the case of an 'impulsive' injection the radiative losses induce a cutoff, and not a break, 
above 7b r - 

EDITOR: PLACE FIGURE |TJ HERE. 

In Fig. 1 we show the results of calculations when the same B = 0.12/iG and 
A^inj = 3 Gyr are assumed for 3 different power-law indices of the electron injection 
spectrum 

Q{l) 0C7- ain jexp(-7/ 7c ), (3) 

where 7 C defines the characteristic maximum energy of the accelerated particles. The 
synchrotron fluxes are normalized to the radio flux from the Coma cluster observed at 
400 MHz. For the assumed magnetic field, the IC 'counterpart' of this normalization point 
in Fig. 1 is found at the photon energy e « 1 MeV, and therefore for different indices 
the IC fluxes predicted in the X-ray region are somewhat different. Because the synchrotron 
flux at a given frequency depends on the overall number of radio electrons and magnetic 
field as L S y n oc iV e i?( 1+Qc ^ 2 , the X-ray fluxes can be increased or decreased assuming slightly 
lower or higher magnetic fields. Thus for the injection spectrum with = 2.2 a field 
somewhat smaller than the B pa 0.12 /iG assumed in Fig. 1 would be needed to explain the 
observed fluxes both in the EUV and in the X-ray regions. In the case of a steep spectrum 
with a in j = 2.6 an increase of the field to B pa 0.15 /iG could reduce the IC flux to within 
the range of the observed 'excess' fluxes in the X-ray region. But, in addition, we have to 
assume a more recent time for the injection of relativistic electrons in the ICM - in order 
to induce the radiative break at higher energies and to explain also the excess EUV flux. 



EDITOR: PLACE FIGURE g HERE. 
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Thus, uncertainties in the reported fluxes and spectral indices of the radiation in 
different energy bands allow some latitude in the predictions of the magnetic field in the 
Coma cluster. In Fig. 2 we show the spectra of IC radiation calculated for two different 
sets of the parameters B, a in j, and At in j. Note that the value of B ps 0.15 /iG needed in the 
case of a steep radio spectrum with a r — 1.3 (i.e. a in j = 2.6) agrees well with the magnetic 
field deduced by Fusco-Femiano et al. (1999) from comparison of the radio fluxes with the 
excess X-ray flux detected by the BeppoSAX detector. A further increase of the magnetic 
field, up to B < 0.2/xG (Rephaeli et al. 1999), would be possible only for a smaller flux of 
the IC X-rays. 

The total energy in relativistic electrons which is needed for the interpretation of the 
excess EUV and X-ray fluxes in Fig. 2 is equal to W c = 6.3 x 10 61 erg and W e = 1.8 x 10 62 erg 
for a in j = 2.3 and a in j = 2.6, respectively. For the angular size of the EUV excess 
FWHM PS 19' x 13' (Bowyer & Berghofer 1998) the mean electron energy density is of order 
1 eV/cm 3 , and exceeds the magnetic energy density wb = 2.5 x 10" 4 (-B/0.1 /zG) 2 eV/cm 3 by 
more than 3 orders of magnitude. For such a high total energy of relativistic electrons, with 
the spectra extending beyond 100 GeV energy (needed for the interpretation of the radio 
fluxes) , the fluxes of IC gamma-rays should extend to the region of high energy 7-rays. In 
the case of power-law spectra in the radio region with a r ~ (1.1 — 1.2), the extension of the 
IC fluxes with the same index from the X-ray to the 7-ray region predicts fluxes observable 
for the future GLAST detector (see solid line in Fig. 2). For steep spectra with a r > 1.3, 
the predicted IC flux would be only marginally observable. However, because of the high 
gas density in the Coma cluster (up to 3 x 10~ 3 cm -3 in its centre, Briel et al. 1992) the 
flux of bremsstrahlung 7-rays should be observed by GLAST in any case (see Fig. 2). 

In conclusion to this section we note that the models invoking a single electron 
population as considered above may encounter a morphological problem. As pointed out 
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by Bowyer & Berghofer (1998) the spatial extent of the EUV emission < 19' is significantly 
smaller than the ~ 30' size of the 30 MHz radio emission ([Henning 1989|) , even though 
the target photon field for IC EUV emission is uniform. This would mean that either the 
ICM magnetic field somehow increases at larger distances from the core of Coma, which is 
very implausible, or that radio emitting electrons are more spread out than less energetic 
EUV emitting electrons. Although the latter possibility cannot be excluded in principle for 
single electron population models, scenarios assuming a different origin of the EUV and 
radio electrons would allow different spatial distributions of these electron populations more 
naturally. In the following section we consider such a model. 



3. IC origin of the EUV flux in high magnetic fields 

The conclusion that the magnetic field in Coma cluster is small, B ~ 0.1 /xG, is 
an unavoidable result of the assumption that the observed excess X-ray flux has an 
inverse Compton origin. However, the values of the magnetic field in Coma deduced from 
Faraday rotation and depolarization measurements are at a level of B > 1 /iG ( [Kim et "aT 



|1990| , [Feretti et al. 199"5| ), which is also typical for other clusters ( |Kronberg 1994| ). 



If this is the case, then an IC origin of the excess hard X-rays in Coma is absolutely 
excluded. Indeed, the target photon field for such IC radiation is well known. Therefore 
the total number of electrons in the energy range 7~ 5 x 10 3 — 10 4 is fixed. Because 

these energies are larger than 7 r .i ~ 3 x 10 3 B~q 5 for fields of the order of 1 fiG (see 
section 2), the synchrotron radiation of those 'IC X-ray' electrons would be at frequencies 
v > 10 MHz, accessible for radio observations. An increase of the magnetic field by one 
order of magnitude compared with what is assumed in Fig. 1 would therefore result in 
overproduction of the observed radio fluxes at least by two orders of magnitude. 
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Thus the IC interpretation of the hard X-ray excess should be abandoned and other 
explanations must be explored, such as bremsstrahlung emission of the suprathermal 
electrons ( [Kaastra et al. 1998| , [Enfilin fc Biermann 1998]) , or perhaps invoking the thermal 
emission of a multi-temperature plasma - with temperatures extending beyond lOkeV. 

However, for the excess EUV flux an IC origin can still not be excluded if in the energy 
region from 7 euv ~ 300 to ~ 3000 the energy distribution of the electrons would decline 
extremely fast in order to reduce the total number of the electrons with 7 > 7 r l by those 
two orders of magnitude relative to a standard power-law ^(7) with a c > 3. Basically, this 
implies a sharp cutoff in ^(7) somewhere in the gap between 7 euv and j T} \. This scenario 
then suggests that the population of relativistic electrons with energies below a few 100 
MeV, responsible for the EUV emission, should be different from the population of electrons 
of higher energies producing the radio emission in Coma. 

Here we suggest that the low energy electrons can be naturally connected with the old 
population of relativistic electrons injected into the ICM more than a few Gyrs ago, and 
perhaps even as early as at the epochs of the initial starburst activity in the Coma galaxies. 

Powerful galactic winds due to early starburst activity in galaxy clusters, which 
appears to be necessitated by the observed high iron abundance, could enrich the ICM 
also with relativistic particles ( |Volk et al. 1996| ) and magnetic fields ( |Volk fc Atoyan 
1999|) . The possibility that low energy electrons responsible for the excess EUV emission 
in a galaxy cluster could be the relics of such a starburst epoch has been considered by 
Sarazin & Lieu (1998). Here we develop a detailed model for a 2-population origin of the 
relativistic electrons in the Coma cluster. The model in particular suggests that the EUV 
emitting electrons have a starburst origin, whereas the radio electron have been produced 
as secondaries during recent epochs. A sharp radiative cutoff - instead of a radiative break 
- below 7 r l will be naturally produced in the spectrum of the relic starburst electrons, 
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provided that the injection/acceleration of these electrons in the ICM has ended at some 
epoch z > 0.05 (corresponding to the last At > 10 9 yr). 



3.1. Model calculations 

Assuming for convenience of calculations an Einstein-de Sitter Universe (Qq = 1), the 
cosmological time is 

t{z) = mtw^ lAx 1010(1 +2r3/2yr - (4) 

chosing H = 50 (km/s)/Mpc for the Hubble constant. Integrating then Eq.(l) over time, 
one easily finds that all electrons which have been injected into the ICM at an epoch z > 
would at present be cooled down, due to IC losses, to energies below 

310 

7cUt ~ (1 + Z )5/2 - 1 • W 

If these electrons were not reaccelerated at some z' < z, their energy distribution should 
have a sharp cutoff at 7 = 7 cu t(z)- Note that for z ~ 1 or larger, 7 cut oc (1 + z)~ h l 2 . This 
dependence readily follows from the comparison of Eqs.(l) and (4), which indicates that 
most of the radiative cooling of the electrons, and thus the formation of the cutoff in the 
energy spectrum, occurs at the early epoch as the injection/acceleration of relativistic 
electrons in the ICM drops. 

The kinetic equation describing evolution of the energy distribution of electrons 
N = iV(7, t) in time reads: 

dN d 

where P = P(7,t) is the energy loss rate and Q = Q(7, t) is the source function of the 
electrons injected into ICM. Assuming that the injection has started at some time to, the 
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general solution to this equation can be written as (|Pacini fc Salvati 1973 , Sarazin 1999 ): 



* 8T 

N(>y,t) = I Q(r 7 ,r)Vdr, (7) 



/to ^ v 7 ' ' ' d'j 

where T 7 = T 7 (t, r) is the energy of an electron 7 at time r < t, which can be found from 
the equation 

^ = -P(T 7I r) (8) 
with the initial condition r 7 (£, r = £) = 7. 

The most important term in P(7, t) is connected with the radiative losses given 
by Eq.(l), first of all with the IC losses of the electrons in the cosmic MBR field, 
w m br oc (1 + zY oc [t/t(0)]~ 8 / 3 (using Eq. 4). Two other terms are connected with the 
interactions of the electrons with the ICM gas. The ionization (Coulomb) losses depend 
very weakly (logarithmically) on the energy 7 of the relativistic electrons, and for our study 
they can be approximated with a good accuracy as[] (|e.g. Ginzburg 1979|) : 



P ion «4.1xl0- 5 n gas yr- 1 . (9) 
The energy loss time for an electron with energy 7 then is: 



7? 



iion = 7/Pon ~ 2.4 x 10 8 — ZT^T 7 yr • (10) 



10~ 4 cm 



-1 



Thus electrons with 7 > 100, which are responsible for the IC EUV emission, will not be 
affected by the Coulomb losses if the ICM gas density n gas (the average during the time 



after injection) would not exceed 10 cm . However at lower energies, depending on n 



gas, 



the Coulomb losses may result in a significant flattening of ^(7) compared with the initial 
injection spectrum. 



1 Note that for a non-ionized medium the ionization loss rate is somewhat smaller, by a 
factor of ~ 2.5. 
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In calculations of the energy distribution of electrons N(^y,t) we take into account also 
the bremsstrahlung losses of the electrons. Note that these losses do not play any noticeable 
role for the formation of iV(7,i) because at all energies 7 they remain significantly smaller 
than the sum of the ionization and IC losses (see Sarazin 1999). However bremsstrahlung 
becomes an important process for production of the 7-ray fluxes (see Sections 2 and 4). 

For the further calculations we approximate the injection spectrum of electrons Q(j,t) 
as a separable function of variables, Q(j,t) — > Q(j) x q(t), with the energy distribution 
(5(7) given by Eq.(3), and with a time profile of the injection rate approximated in the form 

q(t) = 1 , t(z ) <t<t( Zl ) (11) 

[t( Zl )/t] m , t( Zl ) < t < t(0) 

Here the epochs z and Z\ correspond to the beginning and to the end of an efficient 
injection of relativistic particles (accelerated mostly at the early galactic wind termination 
shocks, see Volk et al. 1996) into the ICM, so that the time interval At inj = t(zi) — t(z ) 
has the meaning of a characteristic duration of the electron injection phase ('starburst/wind 
activity'). 

3.2. Origin of the 'EUV electron population 

The electron energy distributions calculated for q(t) declining with different indices 
m are shown in Fig. 3. A stationary injection during At in j = 2 Gyr starting from 
t(zo = 1.5) ~ 3.3 Gyr is assumed. This results in the epoch z\ ~ 0.8 when the injection 
rate has sharply declined. A significant cutoff or a 'knee' in the electron distribution 
above 7 cu t(0-8) = 93 is formed only if the injection rate drops sharply, m > 4. Otherwise 
the electron injection currently remains high enough in order to fill up the knee with 
recently accelerated electrons, resulting in a feature resembling more an ordinary spectral 
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break rather than a cutoff at 7 ~ 7cutO z i)- Then, as it is seen from Fig. 4 (solid curves, 
corresponding to m = 2), in an ICM magnetic field B > 1 /zG the synchrotron fluxes would 
significantly exceed the radio fluxes observed, although the IC flux is still by an order of 
magnitude less than the excess X-ray flux. 

Such an extraordinarily fast decline of q(t) would actually mimic a Gaussian type 
function, with a dispersion a t ~ At- m ^/2 significantly smaller than the time interval 
t(0) — t(zi) between the epoch z > z\ of the maximum injection rate and the present time 
t(0). Note that this kind of injection could be a quite appropriate approximation for the 
electrons injected into the incipient ICM at the epoch of early starburst activity. 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE | HERE. 

The absolute values of the fluxes in Fig. 4 could be changed assuming a different total 
energy for the injected electrons. However the cutoff in the IC radiation would remain 
below the EUV region. Formally, the cutoff in the energy distribution of the electrons could 
be shifted to values 7 cut > 300 needed for an explanation of the EUV flux, if we would 
assume that the phase of active injection of the electrons were continued until very recent 
epochs, 0.05 < Z\ < 0.2, but then it just stopped. However, this kind of assumption could 
hardly find any physical justification, especially if this population of electrons is due to the 
early starburst activity in Coma. 

Much more realistic seems to be a scenario where relativistic electrons driven into 
ICM by the galactic winds at the times of the early starburst activity at z s h > Z\ ~ 1, had 
cooled down to energies 7cut(^i) < 100, remaining however relativistic. The number of these 
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electrons, N e (rel) ~ 10 66 per each 10 61 erg of initial energy, is very large. Then, the passage 
of a strong shock(s) across the cluster, e.g. due to 'explosive' galaxy mergers during the 
formation of the cluster (|Cavaliere et al. 1992| , [Kontorovich 1992|) , would be sufficient to 



reaccelerate these electrons to GeV and higher energies. In the downstream region after the 
passage of the shock the electrons would again be radiatively cooled down. A cutoff in the 
energy spectrum of the electrons will then appear at energies 7 > 500 if such a shock(s) 
occurred at z < 0.2 (i.e. during the last At < 3 Gyr). 

Another interesting possibility to 'bring' the cutoff in the energy distribution of the 
electrons into a 'proper' place is connected with a purely adiabatic energization of that relic 
population of electrons by compression of the ICM gas due to contraction of the cluster. 
An interesting feature of the adiabatic compression is that it does not change the shape of 
the initial energy distribution of the electrons, but just shifts their energy up by a factor of 
a: N^j) = iVo(7 / 'a) / 'a. Therefore, after this kind of energization, the cutoff in the energy 
distribution of the relic population would be conserved, and there is no need to assume that 
the compression has stopped some At > 10 9 yr ago in order to allow the formation of the 
cutoff due to radiative losses (the compression may go on even presently!). For relativistic 
particles the parameter a = (ngas/^1) 1 ^ 3 , where n\ and n gas are the ICM gas densities before 
and after compression. Assuming that the electron injection connected with the starburst 
activity in the cluster was basically over at say z\ ~ 1, for the spectral cutoff formed at 
that epoch we would have 7 cut ~ (60 — 70). Pushing this cutoff up by a factor of a ~ 10, we 
would explain the observed excess EUV flux of Coma by IC radiation of the adiabatically 
reenergized relic population of electrons, but avoiding at the same time an overproduction 
of the radio fluxes at frequencies v > 10 MHz (accessible for the radio observations of Coma) 
in a high ICM magnetic field B > 1 fiG. 

An adiabatic energy amplification factor a ~ 10 implies an increase of the gas density 
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by a factor ~ 10 3 . In principle, a compression to even significantly larger factors is quite 
reasonable in the cooling flows of galaxy clusters like A2199 (see |Kaastra et al. 1999| ). 
Although Coma appears to be a cluster without such a flow, compression of the gas up to 
3 orders of magnitude in its central region could be possible if the formation of the gas 
density profile occurred not earlier than during the last several Gyrs. Note that the angular 
size of the excess EUV emission region ( |Bowyer Berghofer 199"8| ) correlates well with the 
radius of the central core 9 C ~ 10' of the (thermal) X-ray emission region ( [Briel et al. 1992] ), 
where the maximum gas density is as high as ~ 3 x 10~ 3 cm -3 , and the mean gas density in 
the entire EUV emission region can be reasonably estimated as n gas ~ 10 -3 cm~ 3 . Thus, the 
parameter a can be about 10 if the gas density profile in Coma started its formation (due to 
overall contraction, galaxy mergers) from rii ~ 10~ 6 cm~ 3 . This would then imply that the 
ICM gas in the core of the Coma has been increased to the level of 10 -3 cm -3 not earlier 
than during the epochs z < 0.5, because otherwise this initial value of ri\ were smaller than 
the primordial baryonic density. Indeed, given the primordial baryonic density at present 
ribar(O) — 3 x 10~ 7 cm -3 , this density increases as (1 + z) 3 to 10 _6 cm~ 3 at z = 0.5. Thus, 
the compression timescales should then have been less than i(0) — i(0.5) = 6 Gyr. 

EDITOR: PLACE FIGURE | HERE. 

It is important however that an essentially smaller degree of gas compression will 
be needed if we combine this adiabatic compression scenario with the merger shock 
reacceleration at some z' < 1. In Fig. 5 we show an example of this kind of interpretation 
of the excess EUV flux of Coma by IC emission (heavy solid line) of the early population of 
electrons which have been (re) accelerated in ICM by the merger shocks at times z ~ 0.5, 
and later on compressed by a factor a = 3. This would require a compression of the gas 
only by a factor of 27. Note that the gas density at radii 9 ~ 100' (corresponding to 
r ~ 4Mpc) inferred from the (thermal) X-ray fluxes is ~ 10 -5 cm~ 3 (Briel et al. 1992), i.e. 
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the gas compression in the central < 1 Mpc region of Coma up to a factor of ~ 100 is quite 
plausible. 

3.3. Impact of ICM gas 

For calculations in Fig. 5, for the mean ICM gas density averaged in time over the 
last A m 6 Gyr (corresponding to z < 0.5) we have assumed n gas = 10~ 5 cm~ 3 . In that 
case the Coulomb losses of the electrons become noticeable (assuming also that the ICM 
gas was essentially ionized during that time) only at energies 7 < 10, resulting in some 
flattening in the spectra of IC radiation below 0.1 eV (see Fig. 5). Calculations show that 
the Coulomb losses have no a significant impact for interpretation of the EUV flux also for 
n gas ~ 10~ 4 cm~ 3 . However, the assumption that the gas density in Coma was close to the 
current mean n gas ~ 10 -3 cm -3 during the recent epochs z < 0.5 dramatically changes the 
results. Indeed, as follows from Eq.(10), in that case the electrons with 7 < 300 would be 
effectively removed by Coulomb losses. Thus, the IC origin of the excess EUV flux implies 
either that the increase of the gas density in the central region of the Coma has occurred 
on timescales no more than 2 Gyrs (in the framework of the model invoking adiabatic 
energization), or that the relic electrons in the core have been reaccelerated by the shock(s) 
(1-2) Gyrs ago. Note that a shock propagating with a speed of a few 1000 km/s travels 
1 Mpc during a time significantly less than 1 Gyr. In clusters with strong cooling flows, 
where the gas density could reach values well in excess of that in Coma, it would be even 
possible to expect a deficit of the excess EUV emissivity in the centre, if the infall timescale 
there would exceed the Coulomb loss time. 
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3.4. Origin of the 'radio' electrons 

The interpretation of the radio fluxes from the Coma cluster in the case of high 
magnetic fields requires a second component of relativistic electrons injected into the ICM 
at more recent times At < 10 9 yr. Because the overall number of such electrons needed for 
production of the observed radio fluxes is relatively small, the observed excess X-ray flux 
cannot be produced by those electrons (see heavy dashed curve in Fig. 5), and some other 
interpretation (e.g. multi-temperature plasma emission? ) has to be considered. At the 
same time, magnetic fields B > 1 /xG allow a scenario where this second component of the 
electrons in Coma is due to CR protons and nuclei interacting with the gas (via production 
and decay of 7r ± -mesons, pp —> 7r ± —>■ ^ —> e^), whereas such a secondary origin of 
the radio emitting electrons in the case of low magnetic fields B ~ 0.1 /xG would require 
unrealistically high energy in CRs (see |see Blasi fc Colafrancesco 1999|) . The injection rate 



of radio electrons assumed in Fig. 5 could be provided by relativistic protons with a total 
energy M^cr < 3 x 10 62 erg (depending on the CR spectrum at energies E v < 100 GeV). This 
seems to be a quite reasonable value, which can be compared with the energy W e ~ 10 61 erg 
assumed to be injected into Coma in relativistic electrons initially (via galactic winds), 
and with the energy W e ~ 5 x 10 61 erg in the relic electrons presently - i.e. after their 
compression and/or reacceleration. Note that from the point of view of the secondary 
origin of radio emitting electrons in galaxy clusters it may also be meaningful that the 
Coma cluster is characterized by a very high mean gas density up to Mpc spatial scales, 
comparable with the scales of the observed radio emission. The difference in the sizes of 
the EUV and radio halos can be perhaps connected with different spatial distributions of 
high-energy relativistic protons and low-energy relativistic electrons. The reason is that 
the relic population of low energy electrons could be more efficiently compressed (and 
adiabatically reenergized) in the cluster core due to cluster contraction than the secondary 
radio electrons produced currently from very high energy protons. Both these protons and 



-20- 



their secondary electrons can much more effectively propagate across the cluster. 

Thus, for high magnetic field scenario there is no need for current acceleration of 
particles in Coma at all. Moreover, an IC origin of the EUV emission in that case allows 
current acceleration to occur only in some parts of the cluster, but not throughout of its 
volume, because otherwise most of the 'old' population of electrons would be involved in 
acceleration, the cutoff above 7 euv would be removed, resulting in a strong overproduction 
of the radio fluxes. 

4. Gamma-ray emission. 

An inverse Compton origin of the excess EUV emission requires an overall energy in 
relativistic electrons with 7 ~ 300 of order 5 x 10 61 erg in order to produce the observed 
flux. Given the high gas density in Coma, the expected fluxes of the bremsstrahlung 7-rays 
should be significant in both scenarios assuming low or high magnetic fields. The solid 
curve in Fig. 6 shows the fluxes of bremsstrahlung 7-rays produced by the relic population 
of the electrons which explain the EUV flux of Coma in Fig. 5. 

EDITOR: PLACE FIGURE | HERE. 

A remarkable feature of the bremsstrahlung spectrum in Fig. 6 is the very sharp cutoff 
at energies above a few 100 MeV. Comparison of this spectrum with the 7-ray fluxes shown 
in Fig. 2 reveals that the scenario assuming low magnetic fields in Coma, B ~ 0.15 /zG, can 
be easily distinguished from the case of high magnetic field by an essentially more gradual 
decline (if any at all) of 7-ray fluxes above 300 MeV produced by GeV and higher energy 
electrons due to both the bremsstrahlung and the IC processes. An IC origin of the EUV 
emission of Coma predicts in both cases 7-ray fluxes well above the level of sensitivity of 
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the future GLAST detector which thus would help to resolve the problem of low vs. high 
magnetic field in Coma. This possibility to probe the strength of the magnetic field could 
be available also for other clusters with excess EUV emission if their gas density would be 
high. 

Thus, the detection of a sharp cutoff in the spectrum of high energy 7-rays by GLAST 
would confirm the > 1 fiG level of the magnetic field in Coma. We note however that 
an absence of such a cutoff will not yet mean that the magnetic field is low. At energies 
above 100 MeV one can expect also significant fluxes of 7-rays due to the production and 
the decay of 7r°-mesons if the energy content in relativistic protons in the central < 1 Mpc 
region of Coma is high, with Wqr > 10 62 erg. Fig. 6 shows that the total energy flux 
of the bremsstrahlung (solid curve) and 7r°-decay (dotted curve) 7-rays may have even a 
single power-law form if Wcr ~ 3 x 10 62 erg. On the other hand, Fig. 6 shows also that 
in that case the fluxes of 7r°-decay 7-rays should be detectable also by the forthcoming 
imaging atmospheric Cherenkov telescope (IACT) arrays in the region of very high 
energies, E > 100 GeV (hereafter VHE). Note that the upper edge of the range of expected 
sensitivities of such IACT arrays ([see Aharonian et al. 1997 ) shown in Fig. 6 corresponds 



to an extended source with angular size A8 ~ 1°, and the lower edge is for a 'point-like' 
source (with At < 0.1°), and the expected angular size of the Coma cluster in VHE 7-rays 
is < 0.5°. Thus, in the case of a high CR proton content in Coma the fluxes of the 7r°-decay 
radiation, as well as the size of the emission region, could be measured at both high and 
very high energies. It is important to note also that at GeV and lower energies the spectral 
shape of 7r°-decay radiation is not sensitive to the spectral distribution of CRs at very high 
energies. Therefore in principle it should be possible to recover the spectral shape of 7-ray 
fluxes of electronic origin in the (0.1 — 1) GeV region even if VHE emission of Coma were 
not detected. 
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The measurements of VHE 7-ray fluxes can yield a rather accurate value for the mean 
cluster magnetic field, comparing these fluxes with the radio fluxes at high frequencies. 
This will be possible as far as the intensity of the secondary electrons will be then known, 
because they are tightly connected with the 7r°-decay 7-rays. The angular size of Coma 
100 GeV should be comparable with the compact size of radio emission at GHz frequencies, 
with FWHM ~ 15' ( [Kim et al. IMg , |Giovannini et al. 199^ , Peiss et al. Tgg7j) . 
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Fig. 1. — The synchrotron and IC fluxes calculated for the magnetic field in Coma 
cluster B = 0.12 /iG, assuming stationary injection of relativistic electrons during the last 
At in j = 3 x 10 9 yr, and 3 different power-law indices for the source function Q{^) in Eq.(3): 
a- m ] = 2.2 (dot-dashed curves), a-^ = 2.3 (solid), and a in j = 2.6 (dashed). An exponential 
cutoff energy 7 C = 2 x 10 6 is assumed. The source function is chosen so that the radio 
spectrum is normalized to the observed flux at 400 MHz. The compilation of radio data and 
the slope of the hatched region from 30 MHz up to 1.4 GHz, corresponding to a power- law 
index a r = 1.16 for the differential energy flux J(E) = EF(E), are taken from Deiss et 
al. (1997). The hatched regions in the EUV and X-ray domains correspond to the fluxes 
observed by Lieu et al. (1999), and Rephaeli et al. (1999) and Fusco-Femiano et al. (1999), 
respectively. 
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Fig. 2. — The bremsstrahlung and IC radiation fluxes calculated in the case of injection 
of relativistic electrons with a- m j = 2.3 during the last At in j = 3 Gyr assuming B = 0.1 /xG 
(solid curves), and = 2.6, At in j = 1 Gyr assuming B = 0.15 //G (dashed curves). A mean 
gas density n gas = 10" 3 cm" 3 in the ICM is assumed. In the 7-ray region, the expected flux 
sensitivity of the GLAST detector (from Bloom 1996) and the upper flux limit of EGRET 
(Sreekumar et al. 1996) are also shown. 
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Fig. 3. — The energy spectrum of electrons injected at the epoch z > z — 1.5 during 
At in j = 2 Gyr (i.e. z\ = 0.8) calculated for an injection spectrum Q(j) with a in j = 2.1, and 
with a time profile of the injection rate given by Eq.(ll), assuming different decline rates: 
m = 2 (solid), m = 4 (dashed), m = 6 (dot-dashed), and m = oo (dotted curve). In all cases 
the overall injected energy is W e = 10 61 erg. It is assumed also that the gas density n gas in 
the ICM at early epochs z ~ 1 was not much higher than the primordial baryonic density 
at that time, n har (z = 1) ~ 2.4 x l(T 6 cur 3 (adopting n bar (0) ~3x l(T 7 cm~ 3 ). 
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Fig. 4. — The fluxes of synchrotron and IC radiations produced by the electrons shown in 
Fig. 3. A mean magnetic field B = 1 /xG in the cluster is supposed. 
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Fig. 5. — The fluxes of synchrotron (thin lines) and IC (heavy lines) radiations produced by 
the relic population of electrons injected into the ICM at epochs z > 0.5, but reaccelerated 
at z sh ~ 0.5, and later on compressed by a factor a = (n/ni) 1 / 3 = 3 (solid lines). The 
spectrum of IC radiation calculated for the electrons without assumption of later adiabatic 
compression, is shown by the dot-dashed line. The total energy of the electrons explaining 
the EUV flux is W e (inj) = 5.9 x 10 61 erg. For the assumed magnetic field B = 1 fiG, the 
observed radio fluxes should be produced by a second population of electrons (dashed curve), 
e.g. assuming an injection with L in j = 2.9 x 10 43 erg/s, during the recent At = 1.7Gyrs 
(corresponding to z < 0.1). 
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Fig. 6. — 7-ray fluxes expected from the Coma cluster. The solid curve shows the 
bremsstrahlung flux produced by the relic population of electrons which explains the observed 
EUV flux of Coma in Fig.5. Heavy dashed and thin dashed curves show the bremsstrahlung 
and IC fluxes produced by the radio electrons (in the case of B = 1 /xG). Dotted and dot- 
dashed curves correspond to the fluxes of 7r°-decay 7-rays produced by CR protons with total 
energy Wqr = 3 x 10 62 and 3 x 10 61 erg, respectively (for a pr = 2.1). The lower and upper 
edges of the shown range of sensitivities of IACT arrays correspond to a point source and 
an extended (AO = 1°) source, respectively (from Aharonian et al. 1997). For comparison, 
the flux level corresponding to 10% of the Crab Nebula flux (e.g. Konopelko et al. 1999) is 
also shown. 



